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Abstract 


This  is  a  study  of  the  mechanisms  and  results  of  reso¬ 
nance  absorption  in  v/arm  and  cold  plasmas.  Maxwell's  equa¬ 
tions  and  the  plasma  fluid  equations  (neglecting  ion  motion 
and  assuming  wavelike  solutions  in  the  x-direction)  are 
linearized.  The  linearization  is  accomplished  for  a  plasma 
with  a  positive  number  density  gradient  in  the  z-direction. 
Second-order  equations  are  derived  from  the  linearized  set. 
These  second-order  equations  are  differenced  and  solved  in 
the  z-direction  for  the  TM  mode  of  propagation  using  a  two- 
sweep  algorithm  with  zero  and  radiation  boundary  conditions. 
The  characteristics  of  the  field  quantities  are  investigated 
at  various  temperatures.  Further,  (the  plasma  v;ave)  , 
the  temperature  of  the  hot  electrons,  and  the  nonlinear 
Ponderomotive  force  are  calculated.  Finally,  resonance 
absorption  by  the  plasma  is  calculated  and  peak  absorptions 
of  50%  are  observed  for  both  cold  and  low  temperature  warm 
plasmas.  The  validity  of  the  assumptions  is  discussed 
referencing  both  power  and  temperature  concerns.  Possible 
nonlinear  and  time-dependent  modifications  to  the  theory 
are  discussed.  Landau  damping  is  derived  and  its  limita¬ 
tions  considered.  The  two-sweep  algorithm  is  found  to  give 
accurate  results  and  its  amenability  to  computer  application 
makes  it  a  desirable  method. 


I  Introduction 


In  this  thesis  the  propagation  and  absorpton  of  EM 

waves  in  cold  and  warm  inhomogeneous  plasmas  are  studied. 

It  will  be  useful  to  start  with  a  physical  description  of 

the  fundamental  concepts. 

An  electromagnetic  wave  can  freely  propagate  in 

homogeneous  plasma  in  which  the  natural  frequency  of  the 

2  1/2 

plasma,  =  (47rne  /m)  ,  is  less  than  the  frequency  of 

the  wave,  cj.  Here  n  is  the  density  of  the  electrons  in  the 
plasma,  e  their  charge,  and  m  their  mass.  When  the  fre¬ 
quency  of  the  wave  is  less  than  the  plasma  frequency,  the 
wave  does  not  propagate  but  decays  exponentially  in  the 
plasma.  This  exponentially  decaying  wave  is  called  an 
evanescent  wave. 

An  interesting  problem  occurs  in  an  inhomogeneous 
plasma,  then,  when  an  electromagnetic  wave  is  made  to  prop¬ 


agate  from  an  underdense  region  (w  >  co  )  to  an  overdense 

P 

region  (w  <  co^)  .  It  is  the  interaction  of  the  EM  wave  with 


Fig  1.  Density  n^  as  a  Function  of  z  in  the  Case  of 
Linear  Density  Gradient 
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Fig  2.  Geometry  for  the  Incident  Electromagnetic  Wave, 
TM  Mode 


the  plasma  in  this  situation  that  is  the  subject  of  inves¬ 
tigation  in  this  thesis. 

Consider  an  EM  wave  normally  incident  on  a  plasma  in 

which  n  is  increasing  with  z  and  independent  of  x  and  y  as 

shown  in  Figure  1.  For  a  wave  propagating  in  the  z-direction, 

reflection  occurs  at  the  plane  in  the  plasma  vrhere  w  =  cOp. 

This  plane  is  called  the  critical  surface.  For  the  normally 

incident  wave,  the  critical  surface  is  analogous  to  the 

2  2  2 

classical  turning  point  (reflection  point)  where  cos  0  =  Wp/w 
The  mechanism  by  which  energy  is  lost  by  the  electromagnetic 
wave  is  collisional  damping.  This  loss  is  insignificant 
when  the  plasma  temperature  is  above  600  ev. 

Significant  losses  can  occur,  however,  when  the  EM 
wave  is  obliquely  incident  (at  an  angle  3  with  respect  to 
the  normal)  .  This  situation  is  shov/n  in  Figure  2  for  light 
propagating  in  the  TM  mode.  In  this  mode  of  propagation, 
the  EM  wave  has  a  component  of  its  electric  field  pointed 
along  the  plasma  density  gradient.  For  this  reason,  the  EM 
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Fig  3.  Geometry  for  the  Incident  Electromagnetic  Wave/ 
TE  Mode 


wave  interacts  strongly  with  the  plasma  at  the  critical  sur¬ 
face.  This  interaction  leads  to  an  energy  loss  from  the  EM 
wave  that  is  called  resonance  absorption.  A  model  of  this 
process  is  the  specific  objective  of  the  following  work. 

Figure  3  depicts  the  alternative  obliquely  incident  EM 
wave,  which  is  propagating  in  the  TE  mode,  and  is  not  reso¬ 
nantly  absorbed  because  no  component  of  its  electric  field 
lies  along  the  density  gradient.  As  will  be  shown,  this 
wave  is  reflected  with  no  absorption  at  the  classical  turn¬ 
ing  point  which  is  on  the  underdense  side  of  the  critical 
surface.  The  TM  mode  also  undergoes  reflection  at  this 
plane;  however,  the  z-component  of  the  electric  field 
becomes  an  evanescent  wave  which  "tunnels"  into  the  critical 
surface.  Here,  the  EM  wave  becomes  very  large  and  excites 
an  electron  plasma  wave.  The  enhanced  EM  wave  loses  energy 
through  collisions  and  the  plasma  wave  is  damped  without 
collisions  by  interaction  with  electrons  having  a  thermal 
velocity  equal  to  its  phase  velocity.  These  losses  are  the 
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source  of  resonance  absorption.  The  results  of  this  study 
that  demonstrate  these  ideas  are  included  in  Chapter  IV 
along  with  additional  process  descriptions.  The  absorption 
of  the  EM  energy  by  the  electrons  at  the  critical  surface 
and  the  subsequent  transfer  of  this  energy  into  the  random 
thermal  motion  (temperature)  of  other  electrons  will  also  be 
developed  and  discussed  in  Chapter  IV. 

The  linear  variation  of  density  of  a  plasma  models  the 
plasma  surrounding  a  laser  fusion  target  which  has  been 
slightly  heated  to  create  a  thin  enclosing  plasma.  The 
electromagnetic  wave  considered  is  then  the  laser  beam  used 
to  heat  the  plasma.  Resonance  absorption  is  expected  to 
contribute  significantly  to  the  heating  of  the  pellet 
because  collisional  losses  are  quite  ineffective  as  the 
plasma  (target)  temperature  increases. 

The  aim  here  is  to  develop  a  simple  model  of  the 
absorption  mechanisms  and  effects  using  a  minimum  of  com¬ 
puter  space  and  time.  The  effort  will  be  described  as 
follows.  In  Chapter  II  the  equations  to  be  solved  which 
describe  resonance  absorption  in  a  linear  plasma  are  derived, 
as  is  the  term  that  determines  the  damping  of  the  plasma 
waves  (Landau  damping) .  In  Chapter  III  the  numerical  meth¬ 
ods  -for  solving  these  equations  are  discussed,  and  the 
results  are  presented  in  Chapter  IV.  CGS  (Gaussian)  units 
will  be  used  throughout  this  paper  to  enable  comparison  to 
related  research. 
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II  Linear  Modeling 


Two  types  of  plasmas  were  analyzed.  First,  a  cold 
electron  plasma  with  no  electron  pressure  and  then  a  plasma 
of  warm  electrons  with  a  pressure  gradient.  In  each  case, 
the  ions  were  assumed  to  be  immobile  and  an  electron- ion 
collision  probability  was  included.  The  plasma  density 
gradient  necessary  for  resonance  absorption  was  chosen  as 
a  linear  ramp  in  the  z-direction.  The  linear  analysis  was 
achieved  by  using  Maxwell's  equations  and  the  plasma  fluid 
equation  for  electron  momentum.  This  chapter  will  outline 
the  equations  used  and  assumptions  made  in  order  to  obtain 
the  wave  equations  solved  in  the  plasma.  Complicated 
derivations  of  equations  included  in  this  work  are  contained 
in  Appendix  A. 


Linearizing  the  Equations 

Maxwell's  equations  and  the  continuity  and  momentum 
equations  for  electrons  are  (in  Gaussian  units) : 

^  +  V  .  nV  =  0 

dt  — 

dV 

mn(^  +  vV  +  V  •  VV)  +  Vp  =  enE 


1  dB 

^  ^  E  =  -  -  3^ 

477  .  1  dE 

^  ^  i  ^  E  df 


(2.1) 


E  =  477en 


B  = 


5 


further 


j  =  neV 

(2.1) 

P  =  YTn 

where 

V  is  the  electron-ion  (inverse  bremsstrahlung)  collis¬ 
ion  frequency 

n  is  the  electron  number  density 

V  is  the  electron  velocity 
m  is  the  electron  mass 

y  is  the  ratio  of  specific  heats  (3  in  this  one-dimen¬ 
sional  problem) 

T  is  electron  temperature  in  ev 
j  is  electron  current 
B  is  magnetic  field  intensity 
E  is  electric  field  intensity 
e  is  electron  charge 
P  is  electron  pressure 

Linearizing  approximations  were  made  assuming  wavelike 


solutions 

as  follows; 

n  = 

.  -iu)t  + 

no  +  n^e 

iksx 

V  = 

,,  -iojt  +  iksx 
Vi® 

E  = 

So  * 

iksx 

B  = 

5o  ^  ^ 

iksx 

where 


s  is  the  sine  of  the  incident  angle  6 


Hence,  steady-state  solutions  were  the  object  of  the  inves¬ 
tigation  and  oscillations  around  equilibrium  values  the 
result. 


TE  Mode  Solution 

The  transverse  electric  mode  does  not  exhibit  resonance 
absorption  because  its  electric  field  is  perpendicular  to 
the  plasma  density  gradient.  This  can  be  seen  by  examining 
the  wave  equation  for  E^.  To  do  this,  the  above  quantities 
(Eqs  2.2)  were  substituted  into  Eqs  (2.1)  and  resulted  in 
the  following  set  of  equations  for  the  TE  mode. 


E-  =  i)cB^ 


iksE  +  E'  = 

y  y 

-iksB^  +  B^ 


4Trnj^e 
=  -ikeE 


(2.3) 


an^n^  -  anj|^  =  0 

iks (n^  -  n^)  =  0 


where 


a  = 


i47reYT 


cm  (oj  +  iv  )  ' 


'  denotes  differentiation  with  respect  to  z 


and 


e  a  1  -  (1  +  iv/w) 


(2.4) 


Also,  w  is  the  local  plasma  frequency  which  is  given  by: 
^  /  2 

2 

U)  =  - 

p  m 
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From  these  equations  was  derived  the  fact  that  =  0,  and 
that  the  second-order  equation  for  E  is: 

+  k^(e  -  s^)Ey  =  0  (2.5) 

Thus,  there  is  no  plasma  wave  generated  and  no  resonance 
absorption  occurs  in  the  TE  mode.  This  fact  and  Eq  (2.5) 
are  presented  in  Reference  1. 

TM  Mode  Solution 

A  resonance  is  expected  for  the  TM  case  where  the 
electric  field  has  a  component  parallel  to  the  plasma 
density  gradient.  This,  too,  can  be  gleaned  from  an  exami¬ 
nation  of  the  second-order  equations  derived  from  Eqs  (2.1). 
First,  the  equations  were  solved  for  the  TM  mode  v;ith  a 
plasma  number  density  gradient  again  in  the  z-direction 
and  an  electron  temperature.  The  resulting  four  equations 
model  linear  resonance  absorption  phenomena  and  formed  the 
basis  for  this  work. 

ikB  =  E’  -  iksE 
X  z 

-ikeE  =  aiksn,  -  B' 

X  X 

(2.6) 

-ikeEj,  +  a3nj^  -  an|  =  iksB 

iksE  +  E'  =  4TTen, 

X  z  1 

where  6  =  —  n'  and  all  other  variables  have  been  previously 
Uq  u 

defined. 


8 


Consequences  of  the  Density  Profile .  As  mentioned 

earlier,  the  plasma  density  variation  was  chosen  as  a  linear 

ramp  profile.  The  number  density  increased  from  virtually 

19  -3 

zero  at  the  front  face  to  1.12  x  10  cm  two  free-space 
wavelengths  into  the  plasma.  At  that  point,  the  plasma 
frequency  (uj^)  equaled  the  frequency  of  the  incident  .001  cm 
(10  ym)  light.  This  is  the  point  at  which  resonance  (for¬ 
mation  of  plasma  waves)  occurs  and  is  called  the  critical 
surface.  The  linear  variation  in  number  density  made  it 
possible  to  write: 

no  =  n^  (z/L)  (2.7) 

0  =  1/z 

(2.8) 

e  =  1  -  z/(L  +  izv  /w) 

where 

n  is  the  value  of  n  at  the  critical  surface 
c 

V  is  the  value  of  the  electron-ion  collision  frequency 
^  at  the  critical  surface 

L  is  the  distance  to  the  critical  surface 

z  is  the  position  in  the  plasma  (measured  from  the 

front  face) 

Cold  Plasma .  For  the  cold  plasma  case,  the  terms 
involving  T  in  Eqs  (2.6)  are  set  to  zero.  This  is  the  same 
as  ignoring  the  term  "VP"  in  the  momentum  equation  and 
results  in  the  following  equations; 


9 


ikB  =  E'  -  iksE 
X  z 

ikeE  =  B' 

(2.9) 

-ikeE  =  iksB 
z 

iksE  +  E'  =  4Tren, 

X  z  1 

where  all  the  terms  have  been  previously  defined. 

It  is  inconsistent  to  allow  the  electrons  to  have  no 
temperature  and  assume  a  finite  collision  frequency  since 
the  two  are  related  by  (Ref  2:30,  Eq  3.4): 

V/U)  =  24.7  n^/nQ(T^/^)  ;  T  in  ev  (2.10) 

However,  the  v/o)  damping  must  be  included  in  the  equations 
to  avoid  an  infinite  resonance  condition  at  the  critical 
surface.  Later,  v/u  and  Eq  (2.10)  will  be  used  to  obtain 
an  effective  temperature.  Unless  otherwise  stated,  v/w 
was  chosen  small  enough  so  that  only  resonance  absorption 
was  observed, 

Second-Order  Equations 

The  purpose  of  the  preceding  derivations  was  to  obtain 
linearized  equations  which  could  be  easily  converted  into 
wave  equations  for  the  quantities  of  interest.  The  alge¬ 
braic  manipulations  resulting  in  these  equations  are 
included  in  Appendix  A  and  here,  only  the  second-order  equa¬ 
tions  will  be  presented. 

Only  one  second-order  equation  was  necessary  for  the 
cold  plasma.  It  was  written  in  terms  of  E^  and  took  the 
form  (Ref  1:572,  Eq  3.1): 


10 


k^(e  - 


s^)E 


[E^Cilne)  ■]  ' 


=  0 


(2.11) 


For  the  warm  plasma,  it  was  necessary  to  write  two 

coupled  second-order  equations.  Two  sets  of  these  were 

derived  from  Eqs  (2.6).  The  first  set  couples  N,  and  E  : 

X  z 

E^  +  (k^e  -  k^s^)Ej^  =  (4TTe  +  aik)Nj^  -  aikBN^^  (2.12) 

and 

aNl'  +  (a'  -  a6)N'  +  (4Treeik  -  a'3  -  ak^s^)N.  =  -ike'E 
XX  X  z 


while  the  next  couples  B  and  E^ 


E"  +  (-3  + 


i3o<.ks'‘ 


(e  +  aiks  ) 


27)  e;  ^ 


=  (-  ^  +  k^s)B  + 


sBB' 


(e  +  aiks  ) 


and 


(2.13) 


B"  - 


(S'  .t  +  (k^e  -  k^s^jB 

(e  +  aiks'^) 


=  E '  iks 
z 


(ca '  +  aQc  -  ae ' ) 
(e  +  aiks^) 


where 

^  _ iTy 

cmu)  ( 1  +  iv  /w) 

Li 

•  The  warm  plasma  approximation  is  achieved  by  retaining 
the  VP  term  in  the  electron  momentum  equation.  This  resulted 
in  the  terms  containing  v  /w  in  Eqs  (2.6)  which  appears  in  a 
of  Eqs  (2.13).  This  term  (v  /w) ,  which  has  thus  far  remained 

Xj 

undefined,  is  an  attempt  to  model  Landau  damping  of  the 
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plasma  wave  and  will  now  be  explained.  Ginzburg's  second- 

order  equations  (Ref  3:397,  Eqs  20.35)  for  E  and  B  are 

z 

similar  to  Eqs  (2.13)  although  they  do  not  include  Landau 
damping . 


Landau  Damping 

That  there  is  Landau  damping  of  plasma  waves  has  been 
derived  from  the  Vlasov  equation: 


dt 


+ 


V 


Vf  + 


m 


df 

dV 


0 


(2.14) 


There  is  linear  and  nonlinear  Landau  damping.  The  latter 
is  associated  with  particle  trapping  when  n^  =  n^  and  is 
not  within  the  limits  of  linear  theory.  The  derivation  of 
Landau  damping  from  Vlasov's  equation  is  discussed  in  detail 
by  Chen  (Ref  4:213-240)  and  Ginzburg  (Ref  3:122-132),  and 
I  will  follow  Chen's  method  closely.  The  physical  interpre¬ 
tation  of  linear  Landau  damping  was  hinted  at  earlier. 
Simply,  the  plasma  wave  is  damped  without  collisions  by 
exchanging  its  energy  with  particles  satisfying  the  rela¬ 
tion  jv-V  |t  <  ^  =  4  .  These  are  the  particles  in  the 

z  jc  z 

distribution  that  have  not  yet  traveled  one-half  wavelength 
with  respect  to  the  wave.  Hence,  initial  conditions,  such 
as  the  assumption  of  a  Maxwellian  velocity  distribution, 
are  important. 

The  form  of  Landau  damping  employed  in  this  wor)c  was 
derived  in  the  follo\^^ing  way.  First,  wavelike  solutions 
to  the  Vlasov  equation  are  assumed  with  first-order  pertur¬ 
bations  on  the  equilibrium  values.  The  solutions  take  the 


form: 
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f 


f(r,v,t)  =  fg(v)  +  fj^(r,v,t) 


(2.15) 


where 


f ,  ae 


i(k  -  (jjt) 


Vlasov's  equation  then  becomes: 

df 

-  iojf,  +  ikV^f,  =  -  E 

1  2  1  m  z  dV 


(2.16) 


giving: 


ieE  df./dV 
2  0  z 


m  0)  -  k„V 
z  z 


(2.17) 


From  Poisson's  equation  ignoring  E 


ikE^  =  -  477en^  =  -  4T^e 


d^V 


(2.18) 


For  a  one-dimensional  problem,  the  dispersion  relationship 
becomes : 

2 

^  °  ^  (2.19) 

So,  for  large  V  and  small  damping,  the  dispersion  relation 
z 

results  from  an  integration  along  the  real  axis  plus  a 

semicircle  resulting  from  the  integration  around  the  pole 

in  the  complex  V  plane  at  V  =  w/k.  Thus,  the  dispersion 
z  z 

relation  can  be  given  by  (Ref  4:216,  Eq  7-56): 


1  = 


2 

(a1 

2 

U) 


-  u/k  dV^ 


< 

/ 


(2.20) 


where  » 

k 

P  is  the  Cauchy  principal  value  (— |-  in  this  case)  . 
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The  imaginary  part  of  this  dispersion  relationship  account¬ 
ing  for  the  wave  damping  is : 


2  2  ^^0 

1  =  iTTOJ 

p  z  dV 

^  z 


for  a  Maxwellian: 


dV. 


-2V 


/iT  V  +  n‘ 


"^z 

exp(-t3^) 


V 


th 


(2.21) 


2  2  2 

Therefore,  since  prefer  not  to  make  the 

substitution  o  =  u) 

P 


V  3/2  -V 

-im(o/(jj  )  =  1.88  (—-^)  exp(TT— =  v  /w  (2.22) 

P  ^th  ^th  ^ 


where 

has  replaced  as  the  phase  velocity  of  the  wave 

The  wave  equations  for  and  were  presented  earlier. 
The  wave  number  in  the  z-direction  (k^)  for  these  waves  was 
determined  as: 

=  k  /e  (z)  /  /YT/mc^  (2.23) 

where 

k  =  w/c 

The  complete  derivation  of  Eq  (2.23),  including  assumptions, 
is  included  in  Appendix  A.  The  wavelength  and  phase  veloc¬ 
ity  of  the  longitudinal  waves  were  determined  directly  from 
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k  .  The  desired  result,  v  /w  was  then  obtained  from  Eq 
(2.22)  as  a  function  of  z  throughout  the  plasma  and  is 
shown  in  Figure  4.  Note  that  v  /o)  is  temperature  independ- 

Li 

1/2 

ent,  since  in  the  ratio  each  quantity  has  a  T  ' 

dependence.  Past  the  critical  surface,  the  wave  is  expo¬ 
nentially  attenuated  by  e  so  Landau  damping  is  not  tediously 
modeled  there. 

Still,  the  expression  for  Landau  damping  Eq  (2.22) 
fails  when  and  when  e  =  .  4  or  about 

halfway  into  the  plasma.  In  that  region,  v  /m  is  large 
which  makes  for  a  bad  approximation,  since  the  pole  in  the 
complex  integration  (Eq  2.20)  lies  far  from  the  real  axis. 
This,  coupled  with  the  fact  that  the  plasma  wave  is  gen¬ 
erated  closer  to  the  front  boundary  for  increasing  tempera¬ 
tures,  causes  a  breakdown  of  the  warm  plasma  model  at  high 
temperatures.  This  will  be  discussed  in  Chapter  IV. 

Simply,  for  the  bulk  of  the  results  presented  in  this  work, 
v_/u)  was  taken  as  a  constant  in  the  plasma.  Still,  as  will 
be  shown  later,  the  correct  form  of  v  /w  is  imoortant. 

Li 

Having  presented  here  the  equations,  assumptions,  and  some 
background,  the  next  chapter  will  deal  with  the  numerical 
methods  used  to  solve  the  problem. 
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In 
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l.SO  2.00  2. SO  3.00  3. SO 

DEPTH  (CM)  nio'® 


Ill  Numerical  Methods 


Starting  with  the  appropriate  set  of  Eqs  (2.6)  or  (2.9), 
second-order  equations  were  derived  that  described  the 
fields  inside  the  plasma.  For  the  cold  case,  it  was  found 
that  one  equation  for  E^.  could  be  written.  But,  for  the 
warm  plasma,  a  set  of  two  coupled  equations  was  required. 

The  equations,  (2.11)  and  (2.13),  were  solved  numerically 
inside  the  plasma  subject  to  appropriate  boundary  condi¬ 
tions.  The  values  thus  generated  were  used  to  solve  for 
the  other  field  quantities,  N, ,  V  and  the  Ponderomotive 
force.  This  entire  procedure  will  be  described  in  detail. 

Gridding 

The  numerical  solution  began  with  a  choice  of  a  grid- 
ding  scheme.  Due  to  the  resonance  expected  at  the  critical 
surface,  very  fine  numerical  resolution  was  needed  in  that 
region.  However,  the  same  resolution  was  not  necessary  in 
other  parts  of  the  plasma  so  a  differential  gridding 
scheme  was  adopted.  This  affected  the  form  of  the  central 
difference  approximations  to  the  first  and  second  deriva¬ 
tives  in  Eqs.  (2.11)  and  (2.13). 

The  problem  of  derivatives  across  these  nonuniform 
spaces  was  handled  by  adopting  a  technique  used  by  Smith 
(Ref  5:139).  Beginning  with  adjacent  nonequal  mesh  spaces: 

h  <})h 

f - 1 - 1 

n-1  n  n+1 


The  Taylor  series  expansion  of  a  function  A (2)  about  a  point 


"n"  can  be  written  in  two  ways. 

A  . ,  =  A  +  (fhA'  +  4h^4>^A”  +  (terms  of  order  h^) 

n+1  n  2 

(3.1) 

A  =  A  -  hA'  +  -  (terms  of  order  h^) 

n— 1  n  2 

eliminating  A"  gives: 


A' 


n+1 


<t>^A  ,  -  (1  - 

n-1  n 

<}>h(l  +  4>) 


(3.2) 


and  eliminating  A'  results  in: 


A"  = 


^  ,  +  ({)A  .  -  A  —  <j>A 

n+1  ^  n-1 _ n  n 

A^/2  +  <t>h^/2 


(3.3) 


The  error  in  the  second-dorivative  approximation  is  of 

order  h  and  is  the  driving  source  of  numerical  error. 

The  majority  of  the  plasma  was  gridded  with  equal  mesh 

spaces  in  which  case,  Eqs  (3.2)  and  (3.3)  reduce  to  standard 

centered  differences.  Two  mesh  sizes  were  specified:  a 

“  2 

very  fine  mesh  spacing  (on  the  order  of  —  x  10  )  near  the 

critical  surface  to  resolve  the  resonance  peak  and  a  coarser 
gridding  (on  the  order  of  in  the  rest  of  the 

plasma.  The  size  of  the  mesh  spaces  (coarse  and  fine)  was 
determined  somewhat  by  *  ial  and  error.  Each  was  reduced 
until  the  field  values  and  absorption  did  not  change.  The 
coarser  gridding  was  required  to  be  smaller  for  the  warm 
plasma  case  due  to  the  plasma  waves  present  before  the 
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critical  surface.  In  the  region  between  the  two  extremes, 
each  mesh  space  was  doubled  until  the  spacing  was  slightly 
less  than  or  equal  to  the  predetermined  coarse  grid  maximum. 
The  remainder  of  the  grid  was  then  completed  with  equal 
mesh  spaces  of  that  size. 

Differencing  the  Equations 

To  implement  the  sweep  method  employed,  the  value  of 
each  variable  at  mesh  point  "n"  was  assumed  to  be  repre¬ 
sentable  by: 

in  the  cold  case,  and  for  the  coupled  case 


z (n+1) 

‘^^z(n) 

+  LB,  ,  +  M 
(n) 

(n+1) 

=  ”^z(n) 

+  OB,  ,  +  R 
(n) 

where 

K,  P,  Q,  L,  M,  N,  0,  and  R  are  constants 

These  equations  were  used  in  conjunction  with  Eqs  (3.2) 
and  (3.3)  to  difference  Eqs  (2.11)  and  (2.13).  The  objec¬ 
tive  was  to  solve  the  differenced  equations  for  the  values 
of  the  constants  at  each  mesh  point  in  the  plasma  subject 
to  the  correct  boundary  conditions.  Details  of  the  deri¬ 
vation  are  included  in  Appendix  A. 


19 


Application  of  Boundary  Conditions 

The  equations  were  solved  by  employing  boundary  condi¬ 
tions  at  the  front  face  and  past  the  critical  surface  in 
the  plasma.  The  two-sweep  method  began  by  setting  the 
fields  equal  to  zero  past  the  critical  surface  in  the 
plasma.  Two  wavelengths  past  the  critical  surface  was 
chosen  as  the  point  where  the  fields  vanished.  This  dis¬ 
tance  proved  to  be  adequate  to  assure  sufficient  decay 
since  the  wave  was  decaying  exponentially  in  this  region. 
After  applying  this  boundary  condition,  the  constants  in 
Eqs  (3.4)  and  (3.5)  were  evaluated  at  each  mesh  point.  The 
sweep  began  at  the  point  where  the  fields  vanished  and  pro¬ 
gressed  point  by  point  to  the  front  face,  where  the  radia¬ 
tion  boundary  condition  was  applied.  The  radiation  boundary 
condition  forced  conservation  of  the  fields  at  that  point 
and  determined  their  values.  Then  it  was  a  simple  matter  to 
solve  for  the  field  quantities  point  by  point  using  the 
previously  determined  constants. 

Zero  Boundary  Condition.  The  zero  boundary  condition 
took  the  form; 


E  (4X)  =  B(4X)  =  0  (3.6) 

z 

which  was  interpreted  as:  P=Q=L=M=N=O=R=K=0 

Due  to  this  boundary  condition,  the  constants  not  multiply¬ 
ing  a  variable  (M  and  R  from  Eqs  3.4  and  3.5)  were  found  to 
be  identically  zero  throughout  the  plasma.  Therefore,  they 
were  not  necessary  and  the  equations  became: 


^z (n+1) 

=  "^z(n) 

for  the  cold 

plasma,  and 

^z(n+l) 

=  ke  ,  .  + 

2  (n) 

LB,  , 
(n) 

®(n+l) 

=  NE  ,  ,  + 
2  (n) 

OB,  . 
(n) 

(3.7) 


(3.8) 


for  the  warm  plasma.  In  other  words,  the  fields  were 
assumed  to  vanish  two  free  space  wavelengths  past  the  crit¬ 
ical  surface.  This  was  accomplished  by  setting  the  con¬ 
stants  equal  to  zero  at  that  point.  This  made  some  of  the 
constants  unnecessary.  The  remainder  of  the  constants 
were  then  determined  at  each  mesh  point  marching  toward  the 
front  face  of  the  plasma. 

Radiation  Boundary  Condition.  At  the  front  face,  the 
radiation  boundary  condition  was  applied.  From  the  fact 
that : 


E.  =  E-  e 
iric  0  . 

me 


B.  =  e 

me  0.  ^ 

me 


ik(sx  +  cos9z) 


ik (sx  +  cos9z) 


E  e 

sc 


ik(sx  -  cos9z) 


B  =  B„  e 
sc  0 

sc 


ik(sx  -  cos9z) 


(3.9) 


it  can  be  shown  that 


and 


Esc  +  ikeoseE^^  =  0 


B*  +  ikcos9B  =  0 

sc  sc 


(3.10) 
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and  since  E  =  E.  „  +  E„„  and  B  =  B.  +  B 

me  sc  me  sc 

E'  +  ikeosOE  =  2ikcos0E. 
z  me 

(3.11) 

B*  +  ikeosOB  =  2ikcos0B. 

me 

These  (Eqs  3.11)  are  the  radiation  boundary  conditions 
appropriate  at  the  front  face.  They  were  applied  to  the 
electric  field  in  the  cold  plasma  case  and  to  the  electric 
and  magnetic  fields  when  they  were  coupled  in  the  warm 
plasma  case. 

The  radiation  boundary  condition  was  applied  by  solv¬ 
ing  the  Taylor  series  expansion  of  the  field  at  the  front 
face  simultaneously  with  the  wave  equation  Eqs  (2.11)  or 
(2.13).  The  equations  are  each  solved  for  the  first  deriv¬ 
ative  of  the  desired  quantity  (E^  or  B) .  Since  each  is 
treated  similarly,  only  E  will  be  developed  in  this  dis- 
cussion.  The  solution  is  then  substituted  into  the  radia¬ 
tion  boundary  condition  and  that  equation  is  solved  for  E^. 

When  E  has  been  found,  the  forward  sweep  begins  using 
z 

Eqs  (3.7)  and  (3.8).  The  complete  derivation  of  the  appli¬ 
cation  of  the  radiation  boundary  condition  is  contained  in 
Appendix  A. 

The  forward  sweep  completes  the  calculation  of  the 
differenced  variable (s).  Other  quantities  of  interest  were 
then  calculated  using  Eqs  (2.3)  or  (2.6).  Plots  of  these 
quantities  are  contained  in  Appendices  B  and  C  for  both 
cases  and  a  discussion  of  the  results  follows  here. 
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IV  Results 


A  wide  variety  of  conditions  were  investigated  and  a 
large  body  of  results  obtained  from  the  two  computer  codes 
written  to  implement  the  two-sweep  algorithm  just  described. 
In  this  chapter  the  major  results  will  be  presented.  The 
data  from  which  these  results  were  derived  is  contained  in 
Appendices  B  and  C.  The  figures  in  these  Appendices  indi¬ 
cate  the  spatial  variation  of  important  quantities  in  the 
plasma.  The  quantities  investigated  were  the  x  and  z  com¬ 
ponents  of  the  electric  field,  the  magnetic  field,  the 
oscillation  velocity  of  the  hot  (n^)  electrons,  and  the 
Ponderomotive  force  (time-averaged  Lorentz  force  or  radia¬ 
tion  pressure) .  The  collisional  and  collisionless  (Landau 
damping)  absorption  of  the  electromagnetic  waves  energy  was 
also  determined  as  a  function  of  the  angle  of  incidence. 

Some  of  the  results,  such  as  the  percentage  of  the 
laser  energy  absorbed  by  the  plasma  are  independent  of 
power,  barring  instabilities,  while  others  scale  as  the 

electric  field  (n,  and  V  ) .  Ponderomotive  force  is 

1  os 

directly  proportional  to  power  as  is  the  energy  of  the  hot 

electrons  (T„  ev) .  The  results  will  be  discussed  independ- 
n 

ent  of  power  scaling  and  then  the  effect  of  power  will  be 
considered.  Temperature  effects,  such  as  Landau  damping, 
will  be  dealt  with  as  they  apply.  Finally,  the  validity 
of  these  methods  will  be  discussed  referencing  both  power 
and  temperature  concerns. 
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Parameters 


Some  parameters  were  not  varied  throughout  the  analy¬ 
sis.  They  are  now  identified  for  reference  and  comparison 
with  other  works. 


pine  _ 

sin9 

= 

cos6 

X 

B. 

-1 

X  =  .001  cm  (10  pm  free  space  wavelength) 

01  =  1.88  X  10^^  sec”^ 

19  3 

n  =  1.12  X  10  /cm  (critical  number  density) 

G 

k  =  27T/X  =  6283  cm"^ 

L  =  .002  cm  (critical  length) 

It  is  readily  determined  that  P^,  the  incident  power  in 

2 

this  analysis,  was  120  W/cm  .  As  mentioned,  linear  power 
scaling  was  employed  to  determine  values  of  variables  of 
interest  at  current  laser  powers. 


^x  ^ 

The  field  quantities  and  were  calculated  for 
both  cold  and  warm  plasmas  and  are  included  in  the  Appendi¬ 


ces  for  reference.  Although  does  exhibit  a  peak  at  the 
critical  surface,  its  value  there  varies  from  only  2-4  esu 


over  all  the  cases  investigated.  The  transverse  magnetic 


24 


I 


field  was  virtually  unaffected  by  the  resonemce  at  the 
critical  surface,  and  its  peak  in  the  plasma  remained 
approximately  1.5  esu  throughout  the  study.  Since  does 
not  contribute  to  resonance  absorption  other  than  through 
collisions,  and  drifts  involving  the  magnetic  field  were 
not  considered,  these  two  quantities  were  not  of  significant 
interest  and  will  not  be  further  discussed. 

E  (Longitudinal  Electric  Field) 
z 

The  magnitude  and  characteristics  of  E^  are  of  interest 
because  it  is  the  component  of  the  electric  field  that 
excites  the  plasma  wave  and  is  responsible  for  resonant 
absorption. 

Cold  P] asma.  The  longitudinal  field  in  the  cold  plasma 
exhibits  a  narked  resonance  at  precisely  the  critical  sur¬ 
face.  The  magnitude  of  the  field  is  inversely  proportional 

to  the  value  of  v/u).  The  fact  that  ]  E^  ]  av/w  was  derived 

max 

by  Ginzburg  (Ref  3:388).  Some  comparisons  were  made  with 

the  maximum  electric  fields  predicted  by  other  authors. 

For  various  angles  of  incidence,  results  comparable  to 

those  of  Denisov  (Ref  1:574)  v/ere  observed.  For  example, 

for  sin0  =  .3;  v/w  =  .002,  |e  ]  =52  and  for  v/u)  =  .02, 

max 

|e  1  =  5.1  at  the  same  angle  of  incidence.  The  magni- 

max 

tude  of  E  and  its  deocndencies  are  much  different  for  the 
z 

warm  plasma. 

Warm  Plasma .  Linear  Landau  damping,  as  well  as  colli- 
sional  inverse-bremsstrahlung ,  affect  the  longitudinal 
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electric  field  in  the  warm  plasma  regime.  As  noted  earlier, 
temperature  effects  moved  the  resonance  peak  of  toward 
the  front  of  the  plasma.  The  magnitude  of  the  field  at 
resonance  was  also  much  less  than  that  observed  in  the  cold 
plasma.  The  damping  was  more  severe  and  the  resonance 
moved  further  to  the  left  as  temperature  increased.  The 
fact  that  the  field  maximum  decreased  is  not  surprising 
since  the  turning  point  of  the  wave  varies  with  temperature. 
The  wave  turns  (is  reflected)  when: 

e  =  YT/cmoi(l  +  iv  /w)  (4.1) 

Li 

Hence,  for  higher  temperatures  the  wave  is  reflected  closer 
to  the  front  of  the  plasma  (at  a  larger  e)  and  the  evanescent 
region  becomes  larger.  This  allows  less  of  the  evanescent 
wave  to  reach  the  critical  surface  and  initiate  resonant 
absorption.  The  wavelength  of  this  electrostatic  wave  (E^) 
is,  not  surprisingly,  the  same  as  that  for  n^ : 

A^(z)  =  X  /yT/mc^  /  /e (z)  (4.2) 

This  wavelength  (X^)  derived  from  (see  Appendix  A)  was 
observed  in  the  study.  It  follows  from  Eq  (4.2)  that  for 
higher  temperatures,  in  regions  where  e  is  small,  the 
wavelength  is  greater  than  that  of  the  incident  light  and 
X^  clearly  approaches  infinity  at  the  critical  surface. 

For  low  temperature  cases,  the  wavelength  can  be  much 
shorter  than  the  incident  wavelength  very  near  the  critical 
surface.  Examples  of  this  can  be  seen  in  Appendix  C.  For 
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instance,  one  may  compare  the  wavelengths  of  E  near  the 

z 

critical  surface  in  Figures  C-8  and  C-15. 

Spatial  Dispersion.  To  conclude  this  discussion  of 
the  longitudinal  electric  field,  the  effect  of  spatial  dis¬ 
persion  will  be  addressed.  Ginzburg  (Ref  3:389)  suggests 
that  the  use  of  a  local  value  of  e  will  not  be  valid  if  the 
electron  encounters  vastly  different  electric  fields  during 
its  travel  in  the  characteristic  time  under  consideration. 
Ginzburg's  equation  for  the  distance  over  which  decreases 
to  half  its  maximum  value  in  a  cold  plasma  is  (Ref  3:389, 

Eq  20.24) : 

Az  =  L  v^/w  (4.3) 

The  results  obtained  here  for  the  cold  plasma  obey  this 

relationship.  This  accounts  for  gridding  difficulties 

encountered  as  v/cjj  was  reduced  because  in  addition  to  the 

increase  of  jE^l  ,  the  scale  length  given  by  Eq  (4.3) 
max 

decreased.  The  scale  length  is  on  the  order  of  a  Debye 

length  for  the  cold  plasma  when  v/w  =  .002.  The  gradient 

lengths  observed  in  the  warm  plasma  are  clearly  longer. 

Figure  5  shows  for  each  temperature  at  what  power  the  hot 

electrons  (n^)  will  travel  one  scale  length  in  one  period 

(•^)  .  At  powers  higher  than  this,  spatial  dispersion  should 

not  be  ignored.  The  scale  length  referred  to  is  the  same 

as  that  above,  i.e.  the  distance  over  which  E  decreases  to 

z 

one-half  its  maximum  value.  The  figure  shows  both  the  cold 
plasma  case  and  the  warm  plasma  results  for  v  /w  =  .1.  The 

Xj 
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ELECTRON  TENPEROTORE  (KEY) 


Power  at  Which  Spatial  Dispersion  Becomes  Important 
=  Az)  as  a  Function  of  Background  Electron 
Temperature:  Curve  1;  Cold  Plasma  T  =  (24 .7u)/v)  2/3 , 
Curve  2;  Warm  Plasma  v  /oj  =  .1 


results  for  the  warm  plasma  v;hen  v  /u  =  .67  are  not  signif- 

Li 

icantly  different.  Spatial  dispersion  will  not  be  con¬ 
sidered  as  a  limiting  factor  in  the  models  application. 
Still,  it  warrants  more  study,  and  the  above  is  offered  as 
both  a  motivation  and  a  starting  point. 

The  Plasma  Wave  (n^^) 

A  plasma  wave  exists  because  the  electric  field  has  a 
component  (E^)  along  the  density  gradient  in  the  plasma 
causing  charge  imbalance  and,  hence,  a  restoring  force. 

The  electrons  oscillate  at  the  phase  velocity  co/k^,  where 
k  here  is  the  same  as  the  k  derived  earlier  for  E_ .  In 
the  cold  plasma,  the  enhancement  of  n^  at  the  critical  sur¬ 
face  is  so  large  as  to  cause  concern  over  the  linearizing 
approximation  n^  >>  n^.  In  the  warm  plasma,  the  wave  gen¬ 
erally  has  a  smaller  amplitude  and  is  itself  allowed  to  damp 
noncollisionally  by  exchanging  its  energy  with  colder  elec¬ 
trons  in  the  body  of  the  Maxwellian  distribution.  Thus, 
the  wave  is  a  key  performer  in  both  the  warm  and  cold 
scenarios . 

Cold  Plasma .  The  cold  plasma  wave  is  very  much  like 

E  .  In  fact,  no  second-order  equation  was  solved  for  n  , 

Z  X 

and  .it  was  obtained  bv  usina  the  calculated  values  of  E 

-  -  z 

in  an  equation  derived  from  the  four  basic  relationships 
Eqs  (2.9) . 

n,  =  -  c'E  /f4TTe  (4.4) 

X  z 
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Clearly,  n,  behaves  much  like  E  and,  liJie  E  ,  it  is  pre- 
vented  from  becoming  infinite  at  the  critical  surface  by 
collisions.  The  plasma  waves  relationship  to  leads  to 
the  cold  model  predicting  increasing  n^  for  increasing  tem¬ 
peratures  (decreasing  v/tj)  .  Figure  6  shows  the  variation 
of  n^  with  temperature  for  the  cold  plasma,  and  the  warm 
plasma  with  v  /oj  =  .1.  The  figure  demonstrates  this  by 
detailing  for  what  power  and  temperature  the  peak  plasma 

wave  amplitude  (near  the  critical  surface)  equals 
19  -3 

1.12  X  10  cm  (the  value  of  n^  at  the  critical  surface) . 
This  further  serves  to  show  when  the  linearizing  approxi¬ 
mation  n^  >>  n^  is  valid.  The  region  of  validity  is  shown 
on  the  figure  for  both  the  cold  and  warm  plasma  models. 

The  additional  restriction  on  the  validity  of  the  warm 
model  will  Le  explained  later.  The  characteristics  of  the 
warm  plasma  electron  wave  will  now  be  explored. 

Warm  Plasma.  When  temperature  effects  are  introduced, 
is  reduced  greatly  in  magnitude  and  spread  over  a  larger 
portion  of  the  plasma.  For  high  temperatures  large  waves 
exist  far  from  the  critical  surface  and  may  even  impinge 
on  the  front  boundary.  This  translation  of  the  wave  from 
the  critical  surface  to  the  boundary  occurs  because  as  the 
electron  temperature  increases,  the  hot  electrons  created 
at  the  critical  surface  are  allowed  to  transfer  their 
energy  more  readily  to  the  other  electrons  in  the  distri¬ 
bution.  This  suggests  that  a  temperature  dependent  as  well 
as  spatially  dependent  Landau  damping  is  appropriate. 
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POWER  (WFiTTS/CMkk2) 


1 

i 
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III! 


1^1 


GION  0 
VALIDITY  .i;r 
COLD  PLASMVi;2i' 
MODEL  h> 


ELECTRON  TEMPERATURE  (KEV) 

Regions  of  Validity  (Indicated  by  Arrows)  for  the 
Warm  and  Cold  Plasma  Models.  Regions  are  bounded 
by  the  following  curves  representing  limiting 
assumptions:  Curve  1;  cold  plasma  r.]^=no:  Curve  2; 
warm  plasma  (Vj_^/Lo=.l)  n]^=no»  Curve  3;  warm  plasma 


Nevertheless,  as  mentioned  earlier,  the  electron  plasma 
wave  should  be  greatly  attenuated  in  the  vicinity  of  the 
front  face  because  V„  ~  V  there.  The  reason  for  this 

n  U 

(discussed  in  Chapter  II)  is  that  when  V  =  V  ,  the  pole  in 
the  integration  of  Eq  (2-19)  lies  far  from  the  real  axis 
and  contributes  a  large  imaginary  part  (damping)  to  v  . 

This  effect  is  ignored  in  the  derivation  of  Eq  (2.22)  and 
thus  Figure  4  does  not  show  this  large  damping  and  it  was 
not  incorporated  into  the  model. 

In  the  warm  plasma,  as  in  the  cold,  n^  was  derived 
from  the  field  quantities.  In  this  case  E^,  B  and  Eqs  (2.6) 
were  used.  The  equation  derived  was: 

aiks^ 

n,  =  (e;  +  sBV£)/(- -  +  4TTe)  (4.5) 

and  centered  differences  were  used  for  the  derivatives.  In 
addition  to  describing  the  variation  of  the  peak  of  n^  with 
temperature  for  v  /w  =  .1,  Figure  6  shov;s  where  the  velocity 
of  the  hot  electrons  equals  that  of  the  cold.  The  signifi¬ 
cance  of  this  has  been  discussed,  and  the  plot  establishes 
a  lower  bound  on  the  region  of  validity  of  the  warm  plasma 
model.  The  plasma's  absorption  of  the  electromagnetic 
energy  will  now  be  considered. 

Absorption 

Two  methods  of  absorption  of  electromagnetic  energy  by 
the  plasma  were  analyzed.  They  were  collisional  (electron- 
ion  inverse  bremsstrahlung)  and  collisionless  (Landau 
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damping  of  the  plasma  wave) •  Other  absorption  mechanisms 
including  instabilities,  were  ignored. 


The  power  absorbed  was  calculated  by  comparing  the 
reflected  electric  field  at  the  boundary  of  the  plasma  with 
the  incident  field.  The  reflected  fraction  "R"  was  derived 
in  the  following  way.  Since: 

sc  ,  tot  .  inc 

Z  Z  2 


where 

s  c 

E  is  the  scattered  electric  field  at  the  front  face 
z 

is  the  incident  field  (sinG) 
z 

and  E^°^  is  the  value  of  the  field  at  the  front  face 
z 

SC 

The  value  of  E  is  determined  because  the  values  of  the 
z 

total  and  incident  field  are  known.  From  geometry: 


^  /tanB 

X  z 


and  thus 


R 


<S 


> 

z  sc 


>  . 

z  me 


or 


cos 


(4.7) 


(4.8) 


where 

<S  >  is  the  average  value  of  the  Poynting  vector  in 
^  the  z-direction 

This  determines  R,  the  percent  of  the  incident  field 
reflected. 
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The  absorption  of  the  laser  power  varies  with  the  angle 
of  incidence.  As  the  angle  of  incidence  increases,  the 
classical  turning  point  moves  toward  the  front  of  the 
plasma,  and  the  longitudinal  electric  field  travels  further 
through  the  evanescence  region  before  reaching  the  critical 
surface.  For  small  angles,  is  the  dominant  field  and 
electron  motion  is  nearly  parallel  to  the  density  gradient. 
Hence,  the  electrostatic  charge  separation  is  small  (Ref  1: 
574) .  Therefore,  an  angle  of  maximum  absorption  is  achieved 
between  these  two  extremes.  It  is  accepted  theoretically 
(Ref  3)  and  has  been  proven  experimentally  (Ref  6)  that  this 
angle  lies  between  2C°  and  25°  for  the  critical  length  under 
consideration.  Although  the  absorption  of  the  field  in  cold 
and  warm  plasmas  has  somewhat  the  same  magnitude,  the 
methods  of  absorption  are  different.  This  will  now  be 
explained . 

Cold  Plasma.  The  energy  absorption  in  the  cold  plasma 
is  due  solely  to  the  collisional  damping  of  E^.  This  reso¬ 
nance  absorption  is  shown  in  Figure  7  to  have  a  maximum  of 
about  50%  at  18°.  This  is  in  agreement  with  cold  plasma 
work  in  general  and  for  a  low  temperature  warm  plasma 
investigated  by  Forslund  (Ref  7)  using  methods  similar  to 
those  employed  in  this  work. 

Warm  Plasma.  The  warm  plasma  absorbs  energy  from  the 
electrostatic  wave  through  electron-ion  collisions  and 
also  as  the  plasma  wave  is  Landau  damped.  As  can  be  seen 
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from  Figure  8,  the  absorption  due  to  Landau  damping  appears 


to  compensate  for  the  reduction  in  magnitude  of  and  the 
consequential  reduction  in  electron-ion  collisional  absorp¬ 
tion.  Absorptions  of  nearly  50%  were  seen  for  this  low 
temperature  case  (T  =  637  ev)  which  is  similar  to  the  cold 
plasma  results  and  those  of  Forslund  (Ref  7) .  However, 
absorption  results  for  the  warm  plasma  are  highly  sensitive 
to  the  form  of  Landau  damping  used,  especially  at  high  tem¬ 
peratures  . 

Modeling  Landau  damping  as  a  linearly  increasing  func¬ 
tion  of  z  with  a  maximum,  as  it  appears  in  Figure  4,  proved 
unsuccessful.  This  attempt  was  curtailed  when  widely  dif¬ 
fering  absorptions  were  attained  for  only  small  variations 
in  the  slope  of  the  linear  section.  This  can  be  explained 
by  looking  at  the  different  number  of  electrons  (N^)  that 
are  heavily  damped  between  two  cases  (two  different  slopes) 
Then,  absorption  would  be  a  trade-off  between  this  energy 
exchange  with  electrons  and  the  propagation  of  the  wave 
away  from  the  critical  surface.  For  example,  at  large 
angles  of  incidence  and  high  temperatures,  anomalous  absorp 
tion  occurred  due  to  the  closeness  of  the  plasma  wave  to 
the  boundary  and  the  field  structure  there.  In  one  case 
(T  =  5100  ev  and  9  =  1.1  radians)  greater  than  30%  absorp¬ 
tion  was  observed. 

For  the  above  reasons,  it  is  doubtful  that  anything 
concrete  can  be  said  about  absorption  due  to  Landau  damping 
at  high  temperatures  until  some  of  the  assumptions  made  in 
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deriving  it  are  relaxed.  However,  these  r.esults  lead  one 
to  believe  that  Landau  damping  should  be  temperature 
dependent . 


Ponderomotive  Force 

The  Ponderomotive  force  is  a  nonlinear  force  caused  by 
the  gradient  of  the  electric  field.  As  one  can  see  from 
Figure  9,  it  is  large  at  very  low  power  densities  and, 
including  its  effects,  might  be  the  first  thing  one  would 
suggest  as  an  improvement  to  this  steady'-state  model.  The 
force  is  derived  from  the  electron  momentum  equation  by 
Chen  in  Refs  4  and  8.  His  derivation  is  also  included  in 
Appendix  A  of  this  report.  Ponderomotive  force  is  included 
here  as  a  result  derived  from  the  previously  calculated 
field  quantities  and  E^.  The  equation  solved  was: 


NL 


2  2 
u)  VE 

2  16^ 
U) 


(4.9) 


where 

is  the  Ponderomotive  force 
NL 

and  E  =  E^  +  E^  is  2<E^> 

S  X  z 

Thus,  it  acts  on  both  ions  and  electrons,  but  is  mo^'c 

effective  b\'  the  factor  —  ^n  the  el;ctrons. 

me 

The  behavior  of  F.,,  is  similar  in  both  the  cold  and 

NL 

warm  plasma.  In  both,  it  exists  primarily  at  the  critical 
surface.  The  negative  force  in  front  of  the  critical 
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surface  serves  to  pull  electrons  from  the  higher  density 

region  to  the  lower.  This  results  in  the  density  profiles 

observed  in  experiments  (Ref  9)  and  modeled  by  theoreticians 

(Refs  10  and  11) .  Evidently,  the  force  causes  a  depression 

in  the  density  profile  in  front  of  the  critical  surface. 

This  new  profile,  when  allowed  to  smooth,  becomes  a  step 

profile  causing  the  plasma  to  go  from  underdense  (tOp/m  <  1) 

to  overdense  (w  /u  >  1)  in  a  very  short  distance.  This 
P 

effect  is  modeled  successfully  in  Ref  11  for  an  electron 

temperature  of  2  Kev.  They  observed  a  maximum  Ponderomotive 

3  4  2 

force  of  1.25  x  10  ergs/cm  at  120  W/cm  while  for  2.55  Kev 

3  4 

and  V  /w  =  .67,  a  value  of  1.5  x  10  ergs/cm  was  obtained 

J-t 

in  this  work  with  a  linear  profile.  The  Ponderomotive 

force  is  very  large  in  the  cold  plasma  due  to  the  steep 

gradient  in  .  Figure  9  shows  the  Ponderomotive  force  as 

13  2 

a  function  of  temperature  for  a  power  of  10  W/cm  .  The 
plot  also  includes  the  thermal  pressure  force  of  the  plasma 
kTV(nQ) . 

Hot  Electron  Energy  (T^^^) 

The  velocity  of  the  hot  electrons  is  a  concern  because 
they  may  escape  the  plasma  without  depositing  their  energy 
if  they  get  too  hot.  This  would  reduce  enoi'gy  absorption. 

In  fact,  these  suprathcrmal  electrons  have  been  observed 
and  have  been  given  a  great  deal  of  attention  (Refs  12  and 
13)  . 
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In  this  work  the  temperature  of  the  electrons  was  deter 
mined  by  their  oscillation  velocity  (V  ) .  The  V  in  the 
cold  plasma  is  the  familiar  eE/ri(v  -  ie)  ,  while  in  the  warm 
plasma,  additional  terms  appear  due  to  the  pressure  gradient 
For  the  warm  plasma: 


V 

os 


eE 


iksyTn, 


YTnin> 


m(v  -  iw) 


-  iw) 


n^mCvL 


-  iw ) 


yn^t 

-  iw)mnQ 


(4.10) 


where  all  the  terms  have  been  previously  defined. 

The  pressure  related  terms  in  Eq  (4.10)  are  small  correc¬ 
tions  to  the  electric  field  term  and  it  was  discovered  that 
they  could  be  ignored.  Hov/ever,  the  velocities  plotted  in 
Appendices  B  and  C  include  the  effects  of  these  terms,  as 
do  Figures  5  and  6. 

Note  that  Figure  6  shows  at  what  powers  and  tempera¬ 
tures  the  energy  of  the  hot  electrons  is  equal  to  that  of 
the  background  (cold)  electrons.  Since  Landau  damping 
requires  >>  V^j^,  this  curve  provides  an  important  check 
on  the  v;arm  plasma  model. 

Validity  of  Assumptions 

This  is  a  discussion  of  the  major  assumptions  m.ade  in 
this  report.  An  im.portant  assumption  is  that  the  ions  arc 
immobile.  Other  concerns  are  the  validity  of  the  lineariz¬ 
ing  approximations  and  the  effect  of  spatial  dispersion. 
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Both  Ponderomotive  and  Coulomb  forcer,  act  on  the  ions. 

The  Coulomb  forces  are  small  and  can  be  ignored.  The  result 

of  the  Ponderomotive  force  is  the  modification  of  the  plasma 

number  density  profile  mentioned  earlier.  It  has  been 

shown  inter ferometrically  that  the  Ponderomotive  force  can 

modify  the  plasma  density  profile  in  about  a  picosecond 

(Ref  10:184).  This  result  was  obtained  by  illuminating 

glass  microballoons  with  a  Nd-YAG  laser.  Hence,  the  plasma 

21  -3 

was  penetrated  to  a  critical  density  of  10  cm  .  The 

14  2 

experiment  was  conducted  at  a  power  of  10  W/cm  ,  and  the 

profile  modification  was  observed  when  =  F^^.  The 

time  scale  of  interest  in  the  steady-state  problem  developed 

-15 

here  is  on  the  order  of  10  seconds,  somewhat  shorter  than 
that  above.  The  Ponderomotive  force,  although  large,  will 
not  cause  the  ions  tc  move  appreciably  in  10  seconds. 

The  Ponderomotive  force  is  shovm  in  Figure  9  for  a  power  of 

13  2  /i^xmn  J  2  n 

10  W/cm  .  Since  At  =  y  — p —  >>  — p  at  this  power  and  at 

still  higher  powers,  one  can  reasonably  say  that  ion  motion 
is  insignificant  in  one  time  period. 

The  effect  of  spatial  dispersion  was  discussed  earlier 
and  v/ill  be  summarized  here.  It  seems  that  the  cold  plasma 
model  will  definitely  be  affected  by  spatial  dispersion, 
and  the  warm  plasma  model  should  also  incorporate  it  for 
application  of  the  theory  at  powers  higher  than  those  indi¬ 
cated  in  Figure  5.  Still  certain  density  profiles  could  be 
found  for  both  the  cold  and  warm  cases  where  the  effect  of 
spatial  dispersion  would  be  minimized  (lonqer  critical 
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lengths  and  larger  collision  frequencies  in  mediately  come 
to  mind) . 

For  the  purpose  of  this  work,  two  concrete  limits  were 
chosen  on  the  validity  of  the  models.  For  each  temperature, 
these  factors  limit  the  powers  at  which  the  models  can  be 
applied.  The  first  factor  is  the  validity  of  the  lineariz¬ 
ing  approximation  n^^  <<  n^.  For  powers  greater  t.han  that 
at  which  n^  =  n^  this  approximation  is  invalid.  There  is 
no  lower  limit  placed  on  the  powers  for  which  the  cold 
plasma  model  is  applicable.  The  warm  plasma  is  limited  by 
Landau  damping  to  powers  where  the  velocity  of  the  hot 
electrons  is  greater  than  the  velocity  of  the  cold.  This 
is  also  shown  on  Figure  6.  Spatial  dispersion  remains  as  an 


additional  concern. 


V  Conclusions 


The  linear  models  of  cold  and  warm  plasma  resonance 
absorption  developed  here  give  results  comparable  to  those 
obtained  by  other  theoretical  means  and  experiments.  The 
paper  serves  as  a  collection  of  many  facts,  theories,  and 
ideas  about  resonance  absorption.  The  results  indicate 
that  the  numerical  methods  are  sound  and  the  method  was 
found  to  be  conservative  of  computer  time  and  space.  When 
ion  motion  is  ignored,  the  steady-state  method  is  applicable 
to  powers  of  interest  in  today's  fusion  piogram.  Their 
respective  areas  of  applicability  are  indicated  in  Figure  6. 
Correct  modeling  of  Landau  damping  and  spatial  dispersion 
could  strengthen  or  weaken  the  preceding  statements.  Because 
resonance  absorption  is  a  more  efficient  way  to  absorb  laser 
energy  than  simple  collisional  absorption,  lasers  of  fusion 
devices  should  be  positioned  such  that  maximum  resonance 
absorption  is  achieved.  Of  course,  it  is  not  that  simple, 
since  targets  are  spherical  and  uniform  energy  deposition 
is  desired  for  the  implosion.  Also,  the  ideal  angle  will 
change  with  the  density  profile.  Still,  resonance  absorp¬ 
tion  will  play  an  important  part  in  a  laser  fusion  device, 
and  this  work  is  only  a  starting  poinr  in  the  process  of 
applying  the  technique  to  the  machine.  If  these  models 
are  to  contribute  further  in  this  process,  they  might  be 
improved  in  the  following  ways. 
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Current  work  in  the  field  should  give  one  some  idea  as 
to  what  direction  a  follow-on  study  might  take.  In  recent 
works,  there  is  little  mention  of  Landau  damping  while  ion 
motion  and  subsequent  profile-steepening  have  been  investi¬ 
gated.  This  suggests  that  a  nonlinear  analysis  can  be 
used  to  more  exactly  model  the  warm  plasma,  and  that  time 
dependence  is  also  important.  The  large  effect  that  Landau 
damping  had  on  the  absorption  of  electromagnetic  energy  in 
this  study  is  additionally  troubling.  These  suggest  possi¬ 
ble  areas  where  an  improvement  on  the  current  models  can 
be  made. 

An  important  improvement  to  the  warm  model  would  be 
the  ability  to  correctly  model  the  expected  final  distribu¬ 
tion  function.  This  could  be  incorporated  by  using  the 
Vlasov  equation  instead  of  the  momentum  equation.  If  this 
were  done,  the  assumption  of  a  Maxwellian  distribution 
could  be  relaxed  and  a  bump-on-the-tail  distribution  func¬ 
tion  could  be  used.  Then,  Landau  damping  could  be  tackled 
in  a  straightforward  manner  and  some  of  the  problems 
experienced  here  could  be  overcome.  This  would  also  allow 
for  an  assessment  of  the  interaction  between  the  two  tem¬ 
peratures  (relaxation  time)  that  has  not  been  addressed 
here . 

More  simple,  perhaps,  would  be  to  relax  the  assumption 
of  immobile  ions  in  a  step-by-step  manner.  After  one  case 
has  been  studied,  the  initial  profile  can  be  adjusted  by 
applying  the  '  ^motive  force  calculated  in  that 
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steady-state  case  and  moving  the  particles.  The  new  profile 
could  then  be  illuminated  and  the  procedure  repeated. 

A  similar  approach  could  be  used  to  investigate  non¬ 
linear  effects.  First,  the  problem  could  be  solved  by  the 
methods  herein.  Then,  the  values  obtained  could  be  used  to 
solve  Eqs  (2.1)  to  the  next  order.  This  nonlinear  analysis 
might  prove  fruitful  if  more  accurate  results  were  desired 
or  instabilities  were  a  concern. 

It  would  be  useful  to  investigate  the  stability  of  the 
two-sweep  algorithm  employed  in  this  analysis  and  to  deter¬ 
mine  its  applicability  to  related  sets  of  equations. 

In  conclusion,  the  linear  modeling  of  warm  and  cold 
plasma  behavior  gave  results,  including  absorption,  compara¬ 
ble  to  experiments  and  other  worlc  on  resonance  absorption. 
The  behavior  of  the  important  quantities  was  detailed  and 
physically  explained.  Important  consequences  of  resonant 
absorption  wore  shown  and  the  regions  of  theory  validity 
were  established.  Various  ways  to  improve  the  models  were 
suggested . 
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Appendix  A;  Per ivation.s 


Linearized  Equations  from  Fundamental  Equations  TE  Mode 
Equations : 


dv 

+  vV  +  V  •  W)  +  Vp  =  enE 

V  X  E 

= 

1  dB 
c  dt 

V  X  B 

= 

47t  .  1  dE 
c  ^  c  dt 

(A.l) 

V  •  E 

= 

4'nne 

V  •  B 

= 

0 

j  = 

neV 

For  TE  mode 

(see  Figure  3)  assuming  wavelike  solutions 

E  = 

^Oy 

,  „  iksx  -  iu)t 

iy 

B  = 

^  ^  iksx  -  icjt  ^  o  ^  ^  iksx  -  iot 

+  B,  e  +  B-  +  B,  e 

lx  Oz  Iz 

^Ox 

where  quantities  with  subscript  (1)  are  perturbations  about 

equilibrium 

so : 

V  X  E 

= 

-  E  i  +  J-  K  k  ;  where  E  =  E  -  +  E  ,e 

dz  -y  ax  -y  y  -yO  -yl 

and 

1  Hn  1  dO-  •  j.  1  dB^ 

_  i  ^  _  _  1  -Ox  ^  ^  g  i  Oz 

c  dt  c  dt  c  lx  c  dt 

,  u)  „  -iwt 
+  -  B,  c 
c  Iz 
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therefore,  since  V  x  E  =  -  —  , 

c  at 


dz  dz  '"Oy 


i  -3^  +  ik  B,  6'^“'^ 
c  dt  lx 


E  e  Ef. 

dx  y  dx  Oy 


1  ^ 
I  *  i’' 


(A. 2) 


SE  =  B 

y  z 


(A. 3) 


Now  we  must  use 


V  X  B  = 


4tt  ■  1  dE 
c  ^  c  dt 


First,  we  will  find  j  from  the  momentum  equation: 


[en^E^  +  en^^E^e  +  en^^E^e  +  en^E^e 


dY  ^  -imt,  / 
dz  ~  dz  ^"^l® 


this  gives  upon  neglecting  products  of  second-order  pertur 
bations 


dyTn, 
^*^0^0  dT" 


=  0 


,,  -iu)t,  •  ^  .  ,,  -2iut,  , 

mnQVj^e  (v  -  lu))  +  mn^^V^e  (v  -  iw] 


„  -2icjJt  .  „  -iojt  ,  „  -iojt 

=  en,E,  e  +  en,E„e  +  en^E,  e 


i 


h 

» 

1 


neglecting  second-order  terms  and  droppin;.  subscripts  gives 

-ia)t  .  "  _-iu)t  d  _  .-iwti 


=  [en^E^e-^^^  t  en^E^e  ^  YTn,e 

-imt 


•]/m(v  - 


therefore,  since  3 


-ik  +  i4TTe  n. 


cm< 


(u)  +  iv)  /^Y  '''  cm(w  +  iv)dznQ 


i4TreYTdnQnj^ 


YTi4it  e  1 


cm(ai  +  iv )  dz 


This  yields,  when  accounting  for  direction, 


r 

[-iksB^  +  -^p  - 


dng  n^ 


dn. 


ikeE  3  +  ^  dz  n 


-  a 


dz 


-  iks(nj^)i 


this  gives 


dB. 


-  iksB^  +  - 


=  -  ikeE 


(A.  4) 


and  for  equilibrium  values 

dn. 


"1  3  ^ 

^  dz  n^  dz 


0  and  -  iks(nj^)  =  0 


E"  TE  Mode 


TO  find  the  wave  equation  for  Ej,  substituting  Eqs 


(A. 2)  and  (A. 3)  into  (A. 4) 
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io) 


(A. 5) 


-  iks(sE  =  -  ik  E 

y  ik  dz  dz  y  y 

ks^E  -  i  E"  -  keE  =  0 

y  k  y  y 

E"  +  k^ (e  -  s^)E  =  0 

y  y 

Linearized  Warm  Plasma  Equations  from  Fundamental  Equations 
TM  Mode 

To  find  the  four  basic  TM  mode  equations  and  eventually 
second-order  equations,  we  again  begin  with  Eqs  (A.l) 

dV  , 

mn(g^  +  vv  +  V  •  VV)  +  VP  =  enE 

second  order 

neglecting  products  of  other  second-order  terms  and  linear¬ 
izing  for  wavelike  perturbations 

,  dP,.  ,  dP, 

-  n-iojV,  +  vn-V.  +  -  P^ 

0  —1  0—1  m  az  m  dz  ml 


e(no  +  n^) 
m 


<5o  +  Si> 


SO 


5  "o-o 


1  f!o 

prt  dz 


and 


engVi 


2  dP 

—  (n,E„  +  n^Ej  -  -  -  -  iksP, 

m  1—0 _ 0—1  m  dz  m _ 1 

(v  -  iuj) 


from  VxD  = 


iJL  ■  +  i  ^ 

c  ^  c  dt 


(A. 6) 
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(A. 7) 


4tt 


2  dP, 

.  —  (n-E-  +  n  E - ^  P,  ) 

4t  _m _ 1— 0 _ 0—1  edz _ e  1 

c  (v  -  iu)) 


since , 


E  —  Eq  f"  E^e 


i (ksx  -  wt) 


and  E  has  x  and  z  components. 


_  _i^£  iOJ  £ 

c  dt  ~  c  -X  c  -z 


(A. 8) 


where  E  and  E  are  the  perturbed  quantities  (E  and  E  ) 


Since  we  know  that 


giving 


"1^0 


n^yT  dnp 
eng  dz 


(A. 9) 


from  (A. 6),  (A. 7),  (A. 8)  and  (A. 9) 


d(B  +  B  )  . 

iks(Bg  +  B^)k - ^ -  i  = 


(v  -  iw)c;ayenQ  dz 


.  2  /n,YT  dn„ 
47Te  /  1  0 


yn  iksT  ^  „  dn  ^  ^  - 

— -  i  -  ^  -r^  k  +  n,,E^i  +  n.E^k 

e  e  dz  Ox  0  z 


-  ikE  i  -  ikE  k 

X  ^ 


so,  for  the  perturbed  iruantities 
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dz 


-  ikeE  -  aiksn, 
X  1 


-  ikcE^  +  agrij^ 


-  a 


dz 


=  iksB 


(a) 


(b) 


(A. 10) 


now  V  X  E  =  -  —  ^ 
c  dt 


or,  for  the  perturbed  quantities 
dE 


-  iksE^  + 

z  dz 


X 


=  ikB 


(A. Ill 


and  V  •  E  =  4Tre(nQ  +  n^)  yields 


iksE  + 

X 


dE 

-5 —  =  4TTen- 

dz  1 


(A. 12) 


Cold  Plasma:  E^  from  Linearized  Equations 

Now  we  will  derive  the  second-order  equation  for  E^ 
in  a  cold  collisional  plasma  beginning  with  Eqs  (A. 10) 


^  =  ikcE  and 

dz  X 


-  ikcE  =  iksB 


Employing  Eq  (A. 11)  gives 


=  ikB 


differentiating 


-  ikeE„  =  iksB 
z 


yields 


j  dE 

^  _  1  ^  _ z 

dz  ~  s  'dz  dz 
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therefore 


iks^E  -  A  ^ 

z  dz  dz  dz 


■/ike] 


-  ikeE 


multiply  through  by  ik 


3,2(32  _  d  ^  ^  ^ 

z  dz  e  dz  z 


=  0 


clearly , 
E 


_2  ^ 
e  dz 


E^ (£ne) ' 


dl 


=  [E^(£ne)']' 


and  thus 


k^{e  -  s^)E,  +  [E,(Jlne)’]'  +  E"  =  0 

z  z  z 


(A. 13) 


Warm  Plasma:  NI^  and  from  Linearized  Equations 


Equations  for  the  plasma  wave  in  terms  of  the  electro¬ 
static  component  will  now  be  derived  using  Eqs  (A. 10) , 

(A. 11)  and  (A. 12).  Eq  (A. 11)  gives 
sikB  +  siksE 


E'  = 
X 


(A. 14) 


Eq  (A. 12)  gives 


iksE'  +  E"  =  4TTen' 
X  z  1 


(A. 15) 


Eliminating  F/  betv.-cen  Eqs  (A. 14)  and  (A.  15)  and  using  Eq 
(A.lO.b)  for  iksD  gives: 
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1  2 

(-  ikc,E^  +  a'3n^  ~  a  +  iks  E^)ik  +  I.^  =  4Tienj^ 

or 

+  (k^e  -  k^s^)E^  =  (4tre  +  aik)N|  -  aikSN^,  (A. 16) 

Now  for  ; 

differentiating  (A.lO.b)  w.r.t.z 

-  iktE^  -  ike'E^  +  a'Bn  +  aBn'  -  a'n'  -  an"  =  iksB’ 

(A. 16) 

from  (A. 10. a)  and  (A. 12) 

4iTen  -  E' 

B'  =  aiksn  +  £[ - - - -]  (A. 17) 

s 


substituting  B'  from  (A. 17)  into  (A. 16) 

-  ikcE'  -  ike'E  +  a'Bn  +  aBn'  -  a'n'  -  an" 
z  z 


4iTen  -  E' 

=  iks  aiksn  +  c  [ - ] 

s 

and 

aN^  +  (a'  -  aB)N^  +  (4Ticeik  -  a'B  -  ak^s^)N^ 

=  -  ikc'E^  (A. 18) 


Warm  Plasma: 


Now  we 


B  and  E  in 

y  z 


B"  and  E^  from  Linearized  Equations 

shall  derive  second-order  equations  coupling 
the  warm  plasma.  Beginninr;  with  (A.  10.  a) 
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-  ikeE  + 
_ X  dz 

aiks 


(A. 19) 


using  (A. 12)  and  (A. 19)  one  arrives  at; 


ak^s^E  +  aiksE'  = 
X  z 


-  ikcE^  +  dB/dz 


(A. 20) 


where 


differentiating  (A. 20)  w.r.t.z.  and  using  (A. 11)  for 
(ike'  -  ct'k^s^)E  +  (ike  -  ak^s^)  (ikB  +  iksE  ) 


+  a’iksE'  +  aiksE"  =  B" 
z  z 


(A. 21) 


using  (A, 20)  to  eliminate  E  in  (A. 21) 


-  ,  r -  aiksE'  +  B’ ]  - 

(ike'  -  a'k'^s^)  - -  +  (ike  -  ak^s^) 

.  (ike  -  cik  s^)  , 


(ikB  +  iksE  )  +  a'iksE'  +  aiksE" 
z  z  z 


(A. 22) 


Using  (A.  12)  to  eliminate  n  and  from  (A.lO.b)  and  using 

(A. 20)  to  eliminate  E^^  while  employing  (A.  11)  to  eliminate 

E^  from  the  result,  one  obtains: 

aeiks(-  >.iksr')  .  •ir,-).- 

-  ikeL  + - A— y—  + - 

(ike  -  tk‘'s'^)  (ike  -  <k‘^s‘') 


+  aE'  +  ak^sB  +  ak^s^E  -  aE" 
z  z  z 


iksB 
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or 


aE"  + 
z 


=  (iks  -  otk  s)B  - 


ags  B ' 


(e  +  aiks' 


(A. 23) 


Substituting  (A. 23)  into  (A. 22)  for  yields  for  B"  after 
algebra . 

B"  -  (e'  ^  ,1^2^  _  ,^2^2,3 

V  (e  +  aiks^)  / 

=  E^iks  (ca '  +  a$€  -  ac')/(e  +  aiks^)  (A. 24) 


Cold  Plasma ;  Differenced  Equations  Solved  for  the  Constants 
in  the  Plasma 

how  we  will  difference  Eq  (A. 13)  for  E  in  the  cold 

z 

plasma  employing  techniques  described  in  the  paper.  Eq 
(A. 13)  can  be  written  as; 


3 

z  z 


dz 


[E  b  ■ 
z  z  ■ 


+  E" 
z 


0 


where 

,  2  ,  2, 

a^  -  k  (e  -  s  ) 

b^  =  ( f ne )  ' 


thus 


a  E 
z  z 


+  E  c 
z  z 


+  b  E '  +  E " 
z  z  z 


0 


58 


where 


Accounting  for  differential  mesh  spacing  and  substituting 
^n+1  =  PE^  +  q  and  =  Pn-fEn-i  +  qn-1 


'  bp  b  (1  -  (I)'")  p  ,  ^  \ 

n  n  a  a  B 

'  n  n  n  n  n/ 


b  q  q  b  d) 

n^n  _  __n  ,  n^  _  -i-  p 

ot  3  oi  B  ^n-1 

n  nJ  L  n  nJ 


(A. 25) 


where 


Un  =  4)h(l  +  i}*)  and  B^^ 


2  2  2 
(t>  h  ,  d>h^ 
__  + 


therefore, 


- - 0/6 

% 


b  q  q 

n^n  ^n 

a  B 

n  n 


(A. 26) 


(A. 271 


where 


a  +  c  + 
n  n 


b  p 
n^n 


b  (1  - 
n 


When  the  zero  boundary  condition  is  applied  q(N)  is  set  to 
zero,  and  the  remainder  of  the  qj^'s  become  1  by  Eq  (A. 27). 
All  the  Pj^ ' s  are  specified  by  Eq  (A.  26)  and  at  the  front 
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face,  the  radiation  boundary  condition  in  atplied.  How  that 
boundary  condition  is  applied  follows. 

Eliminating  between  the  Taylor  series  expansion  of 
E^  at  the  first  mesh  point,  and  the  wave-equation  (A. 13) 
one  obtains: 


0  = 


2(E2  -  E^l 


-  e  E'  +  a  E,  +  c  E.  +  b  E' 
h  ^  zl  zl 


where  a  ,  b  ,  and  c  have  been  previouslv  defined  and 
z  z  z  - 

means  evaluated  at  the  first  mesh  point.  Using  the  radia¬ 
tion  boundary  condition  for  E'  and  the  fact  that  E2  =  P]^E^  + 
the  above  eouation  becomes 


2  (p  E.  +  q,  -  E,  )  2 

- — - - - - ^(2ikcoseE.  -  ikcosOE.) 

,  2  h  me  1 

h 


+  a  E,  +  c  E.  +  b(2ikcosOE.  -  ikcosOE.) 
z  1  z  1  me:  1 


Solving  for 


^(2ikcos0E.  )  +  b  2ikcosOE. 
h  me  z  me 


2  ,  2ikcos0  ,  ,  -1  n 

-  +  a^  +  c  -  b^ikcosO 

,2  h  z  z  z 


(A. 28) 


Warm  Plasma :  Pi f f crcnce  Equat ions  Solved  for  the  Constants 
in  the  Plasma 

The  first  r; U~p  in  the  solution  process  is  to  solve  for 
the  constants  relatin;)  adjacent  field  quantities.  Eqs  (A.  23) 
and  (A. 24)  were  solved  in  general  as  follows: 
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I  n»'' 


r 


+  b^B*  +  Cj^B  -  d^z'  +  e^z 


(A. 29) 


•"  +  b  z'  +  c  z  =  d2B’  +  e2B  ;  z  =  (A. 30) 


taking 


,  =fB+gz+h^ 

^n+1  n  n  n  n  n 


D  =kB+^z+P^ 

®n+l  n  n  n  n  n 


When 


differential  mesh  spacing  is  applied  to  the  derivatives 


(A. 29)  becomes; 


b,k_  b,  (1  '  i>  ) 


^  "D-  -  A 


/a,  b, i  d^g  d  (1  -  e^)  V 


-  ^  ^  ^  X  ^  /Vl  ^  X 


^  3  b-i 

"  X  ^n  ’’  X  ^n  A 


(A. 31) 


where 


=  ($h(l  +  4>)  and  D  = 


4>^h^  ^ 

'~T~  2 


Similarly 


from  (A. 30 ) ; 


^9  9  *^9  ^9  ^9  \ 

(1  -!})  )-  ^  k  +  —  f  +  ~r  t  -  B 
^  An  Dn  An  2/n 


D 


^n 


D 


^2 

^  X  '^n 


9 

X  ^  "  <^2 


X  i  )2 

A  n/  n 


^2'^  ,  ^2  .2)  , 
“X  X  Tn-l  ^ 


X  Vl 


(  ^2^n  _  ^  r^  ' 

V  D  A  ^n  A  ^ni 


(A. 32) 


Eqs  (A. 31)  and  (A. 32)  can  be  written  as 


QB  +  Rz 
n  n 


s,  B  .,  +  s«z  -,  +  s~, 
1  n-1  2  n-1  3 


and 


(A. 33) 


TB  +  uz  =  V.B  ,  +  V^z  ,  +  respectively, 
n  n  ln-1  2  n-1  3  ^  ■' 


Since  we  knew 


B  =  k  ,B  ,+?,  ,z  1+p  T  and 
n  n-1  n-1  n-1  n-1  '^n-l 


z  =  f  ,B  T+g  1  z  1+h  , 

n  n-1  n-1  n-1  n-1  n-1 


Eqs  (A. 33)  can  be  solved  simultaneously  for  the  constants 
k,  £,  p,  f,  g  and  h. 

QV  -  Ts 


B  =  and  z 

n  Qu  -  RT  n 


Qu  -  TR 


where  s  and  V  are  the  right-hand  sides  of  Eqs  (A. 33) 
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clearly 


n-1 


^n-1 


QV^  - 

Ts 

QU  - 

TR 

QV2  - 

TS2 

Qu  - 

TR 

1  1 

>  ' 
o 

Ts^ 

(A. 34! 


'n-1 


Qu  -  TR 


and 


-  P-Vi 


n-1 


“n-1 


Pn-1 


Qu 

-  TR 
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Two  equals j-t-'i. -3  - 

ticn.  They  are  derived  as  was  r.q  (h.23,  but  they  involve 


both  and  B^. 
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(A. 37) 


where 


is  the  first  mesh  interval 

is  part  of  Z2  (i.e.  Z2  =  fj^B^  +  +  H^) 


Eqs  (A. 36)  and  (A. 37)  can  be  written  as 


W2B1  +  >^2^1  ~  '^2 


Since  all  the  constants  are  known  z^  and  are  obtained 
simultaneously . 


k„  from  E" 
z  _ 

It  is  possible  to  derive  the  wave  number  and  wavelength 
in  the  plasma  from  the  second-order  equation.  The  equation 
for  in  the  wai'm  plasma  is; 
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fk^s 


iks\ 
a  / 
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(A. 23) 


for  small  a  ignoring  all  terms  except  for  Eq  (A. 23) 
becomes  a  simple  differential  equation  with  an  oscillating 
part  determined  by  the  coefficient  of  E^,  which  is: 


k 


2 


+  iv^/w) 
2 

yT/mc 


(A. 38) 


2 

for  V  /i.o  small  and  neglecting  s  in  Eq  (A.  38)  the  wave 

Li 

number  of  the  oscillating  function  is  the  square  root  of 
the  right-hand  side. 


k^^^  =  k  /e  (z)  /  /yT/mc ^  (A. 39) 

The  wavelength  is  clearly  determined  also. 

Ponderomotive  Force  (F  )  (Ref  4:257-58) 

_ _  NL 

Beginning  with  the  electron  equation  of  motion: 

m  ^  -  e[E(r)  +  VxB(r)]  (A. 40) 

assuming  an  electric  field  of  the  form 

E  =  E  (r) coswt 
s 
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neglecting  VxB  of  Eq  (A. 40) 


mdVj^/dt  =  -  eE(rQ) 

=  -  (e/mo))  EgSinwt  =  dr^^/dt 

and  integrating  over  time: 

2 

6r^  =  (e/moj  )  E^coswt 

Since  r^  ir  the  initial  position  of  the  particle. 
Expanding  E(r)  about  r^  one  obtains; 


E(r)  =  E(rQ)  +  (6r^ 

from  Maxwell's  equation, 
VxE  =  -  dB/dt 


V)E 


r=r. 


can  be  found  as 


-  -  VxE 

(jJ  s 


sinwt 


r=r. 


(A. 41) 


(A. 42) 


(A. 43) 


now  the  second-order  part  of  Eq  (A. 40)  is 

mdV^/dt  =  -  e[(6r^  •  V)E  +  V^xB^]  (A. 44) 

Substituting  (A. 41),  (A. 42),  and  (A. 43)  into  (A. 44)  gives 

after  time  averaging 
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mm 


1 

"2  2  ^  ^^'s 


V)E^  +  EgX(7xEg) ] 


=  F 


NL 
(A. 45) 
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1 


reducinn  by  o;'.pandinc;  the  doubJi-  ;  :  .  i^’.ct 


F 


NL 


1  2  „ 
1  ^rv  9 


mio 


which  is  the  force  on  a  single  electron.  The  force 
is  found  by  multiplying  (A.  46)  by  n^.  In  terr-.s  of 
becomes : 
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NL 
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2,  2 

(jJ  /u 
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z/L  in  our  plasma: 
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(A, 46) 

density. 

this 

P 

(A. 47) 


(A. 48) 


Appendix  B:  Cold  Plasir.a  Data  Curves 


This  appendix  contains  quantities  of  interest  as  a 
function  of  position  in  a  cold  plasma  for  a  power  of  120 
watts/cm^.  Power  scaling  relationships  are  given  in  the 
body  of  the  paper.  The  curves  are  grouped  according  to  the 
collision  frequency  used.  Each  group  contains: 

(E  longitudinal)  in  e.s.u. 

E^  (E  transverse)  in  e.s.u. 

the  hot  electron  number  density  (cm  ^) 
the  Ponderomotive  force  per  cubic  centimeter 

Nil 

^transverse  velocity  of  the  hot  electrons  due  to 
E^  in  cm/sec 

V,  •.  the  velocity  of  the  hot  electrons  due  to 

longitudinal  ■' 

E  in  cm/sec 

B  the  magnetic  field  in  Gauss 

Note  that  used  in  the  paper  is  the  square  root  of  the 

sum  of  the  squares  of  and  .  Since 

each  group  of  figures  is  the  same  save  the  value  of  v/o, 
collective  figure  titles  will  be  given  and  the  figures  sim¬ 
ply  numbered.  The  following  are  the  titles: 

Figures  B-1  -  B-7;  Data  for  v/w  =  .002 
Figures  B-8  -  B-14;  Data  for  v/w  =  .008 
Figures  B-15  -  B-21;  Data  for  v/o  =  .02 

Please  note  that  the  legend  of  each  figure  contains  addi¬ 
tional  information. 
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Appendix  C:  Warm  Plasma  Data  Curves 


This  appendix  contains  quantities  of  interest  as  a 

function  of  position  in  a  warm  plasma  at  a  power  of  120 
2 

W/cm  .  Power  scaling  relationships  are  given  in  the  body 
of  the  report.  The  curves  are  grouped  according  to  the 
electron  temperature  and  Landau  damping  coefficient  (L.D.C. 
or  v^/w)  used.  Each  group  contains: 

ij 

(E  longitudinal)  in  e.s.u. 

Ejj  (E  transverse)  in  e.s.u. 

Nj^/Nq  or  the  ratio  of  the  hot  electron  number 

density  to  the  initial  number  density  or  to  the 

.  19 

value  of  Nq  at  the  critical  surface  (1.12  x  10  ) 

F.„  the  Ponderomotive  force  per  cubic  centimeter 
NL 

V.  the  velocity  of  the  hot  electrons  in  the 

transverse  •' 

x-direction  cm/sec 


V,  • j 1  the  velocity  of  the  hot  electrons  in  the 
longitudinal 

z-direction  cm/sec 


B  the  magnetic  field  in  Gauss 


Note  that  used  in  the  report  is  the  square  root  of  the 

sum  of  the  squares  of  V.  and  V,  • j  •  t*  Since 

^  transverse  longitudinal 

each  set  of  figures  is  the  same  except  for  the  value  of 
and  the  temperature,  collective  figure  titles  will  be  given 
and  the  figures  simply  numbered.  The  following  are  the 


titles; 


.0005 


2 

Figures  C-1  -  C-7;  Data  for  Vj^/co  =  .1  and  T/mc  - 

2 

Figures  C-8  -  C-14;  Data  for  Vj^/w  =  .-9  and  T/mc  =  .00125 

2 

Figures  C-15  -  C-21;  Data  for  =  .1  and  T/mc  =  .005 

2 

Figures  C-22  -  C-28;  Data  for  Vj^/to  =  .1  and  T/mc  =  .05 

2 

Figures  C-29  -  C-35;  Data  for  Vj^/co  =  .67  and  T/mc  =  .00125 

2 

Figures  C-36  -  C-42;  Data  for  =  .67  and  T/mc  =  .005 

2 

Figures  C-43  -  C-49;  Data  for  =  .67  and  T/mc  =  .01 

2 

Figures  C-50  -  C-56;  Data  for  Vj^/u)  =  .67  and  T/mc  =  .1 

Please  note  that  the  legend  of  each  figure  contains  addi¬ 
tional  information. 
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